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>v Abstract 

Asteroseismology is the determination of the interior structures of stars by using their 
oscillations as seismic waves. Simple explanations of the astrophysical background and 
some basic theoretical considerations needed in this rapidly evolving field are followed 
by introductions to the most important concepts and methods on the basis of example. 
Previous and potential applications of asteroseismology are reviewed and future trends 
are attempted to be foreseen. 
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Introduction: variable and pulsating stars 

Nearly all the physical processes that determine the structure and evolution of stars occur 
in their (deep) interiors. The production of nuclear energy that powers stars takes place 
in their cores for most of their lifetime. The effects of the physical processes that modify 
the simplest models of stellar evolution, such as mixing and diffusion, also predominantly 
take place in the inside of stars. 
\Q ', The light that we receive from the stars is the main information that astronomers can 

l/-) | use to study the universe. However, the light of the stars is radiated away from their 
surfaces, carrying no memory of its origin in the deep interior. Therefore it would seem 
that there is no way that the analysis of starlight tells us about the physics going on in 
the unobservable stellar interiors. 

However, there are stars that reveal more about themselves than others. Variable stars 
are objects for which one can observe time-dependent light output, on a time scale shorter 
than that of evolutionary changes. There are two major groups of variable star, extrinsic 
and intrinsic variables. 

Extrinsic variables do not change their light output by themselves. For example, the 
light changes of eclipsing binary stars are caused by two stars passing in front of each other, 
so that light coming from one of them is periodically blocked. The individual components 
of eclipsing binary stars are not necessarily variable. By analysing the temporal light 
variations and orbital motion of eclipsing binaries, one can determine their fundamental 
properties, and by assuming that their components are otherwise normal stars, determine 
fundamental properties of all stars, most importantly their masses. In this way, stars and 
stellar systems can be understood better. 

Intrinsic variables, on the other hand, change their light output physically. Supernovae, 
which are stellar implosions/explosions, can become brighter than their host galaxies 
because of the ejection of large amounts of material. Even more revealing are stars that 
vary their sizes and/or shapes: pulsating variables. 
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Table 1: Selected classes of pulsating star 



Name 


Approx. Periods 


Discovery / Definition 


Mira variables 


100 - 1000 d 


Fabricius (1596) 


Semiregular (SR) variables 


20 - 2000 d 


Herschel (1782) 


5 Cephei stars 


1 - 100 d 


1784, Pigott, Goodricke (1786) 


RR Lyrae stars 


0.3 - 3 d 


Fleming (1899) 


S Scuti stars 


0.3 - 6 h 


Campbell & Wright (1900) 


{3 Cephei stars 


2 - 7 h 


Frost (1902) 


ZZ Ceti stars (DAV) 


2 - 20 min 


1964, Landolt (1968) 


GW Virgims stars (DOV) 


5-25 min 


McGraw et al. (1979) 


Rapidly oscillating Ap (roAp) stars 


5-25 min 


1978, Kurtz (1982) 


V777 Herculis stars (DBV) 


5-20 min 


Winget et al. (1982) 


Slowly Pulsating B (SPB) stars 


0.5 - 3 d 


Waelkens & Rufener (1985) 


Solar-like oscillators 


3-15 min 


Kjeldsen et al. (1995) 


V361 Hydrae stars (sdBVr) 


2 - 10 min 


1994, Kilkenny et al. (1997) 


7 Doradus stars 


0.3 - 1.5 d 


1995, Kaye et al. (1999) 


Solar-like giant oscillators 


1 - 18 hr 


Frandsen et al. (2002) 


V1093 Herculis stars (sdBVs) 


1 - 2 hr 


Green et al. (2003) 


Pulsating subdwarf star (sdOV) 


1-2 min 


Woudt et al. (2006) 



The first pulsating star was discovered more than 400 years ago. In 1596 David Fabri- 
cius remarked that the star o Ceti (subsequently named "Mira", the wonderful) disap- 
peared from the visible sky. About 40 years later, it was realized that it did so every ~ 1 1 
months; the first periodic variable star was known (although we know today that in this 
case, the term "periodic" is not correct in a strict sense). 

In 1784 John Goodricke discovered the variability of 5 Cephei, and in 1914 enough 
evidence had been collected that Harlow Shapley was able to demonstrate that the vari- 
ations of 5 Cephei and related stars (also simply called "Cepheids") was due to radial 
pulsation. Also in the teens of the previous century, Henrietta Leavitt pointed out that 
the Cepheids in the Small Magellanic Clouds follow a period-luminosity relation, still one 
of the fundamental methods to determine distances in the visible universe - and one of 
the major astrophysical applications of pulsating stars. 

With the ever increasing precision in photometric and radial velocity measurements, a 
large number of groups of pulsating star is nowadays known. Figure 1 shows theoretical (in 
the sense that the logarithm of the effective temperature is plotted versus the logarithm 
of the stellar luminosity) HR Diagrams containing the regions in which pulsating stars 
were known some 40 years ago and today. Table 1 gives a rough overview of the classes 
of pulsator contained in Fig. 1. 

The different types of pulsator have historically been classified on a phenomenological 
basis. The separation between those types has usually later turned out to have a physical 
reason. The individual classes are different in terms of types of excited pulsation mode 
(or, less physical, pulsation period), mass and evolutionary state, hence temperature and 
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Figure 1: Theoretical HR diagram schematically showing locations of selected confirmed 
types of pulsating star. Parts of model evolutionary tracks for main sequence, horizontal 
branch and post-AGB stars are shown as dashed-dotted lines for orientation. Left: known 
classes of pulsating star about 40 years ago. Right: classes of pulsator known to date. The 
names of the different groups are comprehensively listed in Table 1. Note the enormous 
increase in the number of classes of pulsating star and the improved knowledge about 
their loci in the recent past. 



luminosity. The names of these classes are assigned either after the name of a prototypical 
star or give some description of the type of variability and star. 

It must be pointed out that the present overview does by far not contain all types 
and subgroups of pulsating star that have been suggested. The Cepheids are subdivided 
according to population and evolutionary state, into 5 Cephei, W Vir, RV Tau and BL 
Her stars. Jeffery (2008) proposed a number of types of evolved variable, there are the 
Luminous Blue Variables, and there may be new classes of white dwarf pulsator, oscillating 
red and brown dwarfs, etc. Furthermore, some of the instability domains of different 
pulsators overlap and indeed, some objects called "hybrid" pulsators that show oscillations 
of two distinct types, have been discovered. Also, the instability boundaries of some of 
these variables may need to be (considerably) extended and/or revised in the near future. 
For instance, there may be supergiant SPB stars, and solar-like oscillations are expected 
in all stars having a significant surface convection zone. 

Whereas the classification of and distinction between the different classes of pulsating 
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star, that are historically grown and modified designations, can in some cases be called 
arbitrary today, one recognizes an important fact: pulsating stars populate almost the 
entire HR diagram, and this means that they can be used to learn something about the 
physics of most stars. 

Astrophysical background 
Driving mechanisms 

What can make a star oscillate? After all, stars are in hydrostatic equilibrium: the gravita- 
tional pull on the mass elements of normal stars is balanced by gas pressure. If something 
would hit a star, the inwards moving regions will be heated, and the increased heat loss 
damps the motion. Consequently, self-excited pulsations require a driving mechanism 
that overcomes this damping and results in a periodic oscillation. 

Four major driving mechanisms have been proposed. The e mechanism (Rosseland 
& Randers 1938) assumes a variation in the stellar nuclear reaction rate: if a nuclear 
burning region is compressed, the temperature rises and more energy is produced. This 
gives causes expansion, the pressure drops, and so does the energy generation: the motion 
is reversed and oscillations develop. The e mechanism (where e is the usual designator for 
the nuclear reaction rate in formulae), that operates similar to a Diesel engine, has been 
proposed for several different types of pulsating star, such as our Sun and pulsating white 
dwarfs, but observational proof for oscillations driven by it is still lacking. 

Considerably more successful in explaining stellar oscillations is the n — 7 mechanism 
(Baker & Kippenhahn 1962 and references therein). In layers where the opacity k increases 
and/or the third adiabatic exponent T 3 decreases with increasing temperature, flux coming 
from inner layers can be temporally stored. Such layers in the stellar interior are generally 
associated with regions where (partial) ionization of certain chemical elements occurs. 

The energy accumulated in this layer during compression is additionally released when 
the layer tries to reach its equilibrium state by expanding. Therefore, the star can expand 
beyond its equilibrium radius. When the material recedes, energy is again stored in the 
stellar interior, and the whole cycle repeats: a periodic stellar oscillation evolves. This 
mechanism is also called the Eddington Valve, and it explains the variability of most of 
the known classes of pulsating star. 

The classical pulsators in the instability strip, ranging from the 5 Cephei stars to the 
RR Lyrae stars and the S Scuti stars draw their pulsation power from the Hell ionization 
zone, whereas the oscillations of the roAp stars are believed to be excited in the Hi and 
Hel ionization zones, those of the Mira variables in the Hlionization zone, and those of the 
(3 Cephei and SPB stars are triggered in the ionization zone of the iron-group elements. 

A very similar mechanism, in the sense that it is also due to a region in the star 
behaving like a valve, is convective blocking (or convective driving). In this scheme (e.g., 
Brickhill 1991), the base of a convection zone blocks the flux from the interior for some 
time, releasing the energy stored during compression in the subsequent expansion phase. 
The pulsations of white dwarf stars of spectral types DA and DB as well as 7 Doradus 
stars are thought to be excited (at least partly) via this mechanism, that may also be of 
importance in Cepheids and Mira stars. 
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Finally, the pulsations of the Sun and solar-like stars, that are intrinsically stable and 
therefore not called self-excited, are stochastically excited due to turbulence in their surface 
convection zones. The vigorous convective motion in the outer surface layers generates 
acoustic noise in a broad frequency range, which excites solar-like oscillation modes. Due 
to the large number of convective elements on the surface, the excitation is of random 
nature, and the amplitudes of the oscillations are temporally highly variable. 

Given the physical nature of these driving mechanisms, the existence of the different 
instability domains in the HR diagram (cf. Fig. 1) easily follows. A star must fulfil certain 
physical conditions that it can pulsate, as the driving mechanism must be located in a 
specific part of the star to give rise to observable oscillations. More physically speaking, 
in the case of self-excited pulsations the driving region must be located in a region where 
the thermal, and/or convective time scale closely corresponds the dynamical (pulsational) 
time scale. 

The consequence of the previous requirement is a constraint on the interior structure 
of a pulsating variable: if the instability region of some class of pulsating variable is 
accurately known, models of stars incorporating its excitation mechanism must be able 
to reproduce it. In this way, details of the input physics describing the interior structures 
of stars can be modified to reflect the observations. 

However, this is not the only method available to study stellar structure and evolution 
from stellar pulsation. Models also need to explain the oscillation properties of individual 
stars. We are fortunate to be in the presence of stars having very complex pulsation 
patterns, multiperiodic radial and nonradial oscillators. The research field that determines 
the internal constitution of stars from their pulsations is called asteroseismology. 

It is now due to make clear that the star whose interior structure is best known is the 
star closest to us, the Sun. Its surface can be resolved in two dimensions and millions 
of pulsation modes can be used for seismic analyses. The related research field is helio- 
seismology, and it has been extensively reviewed elsewhere (e.g., Christensen-Dalsgaard 
2002, Gizon et al. 2010). The present article will not touch upon helioseismology. 

Asteroseismology 

The basic idea of asteroseismology (see Gough 1996 for a discussion of why "astero") is 
analogous to the determination of the Earth's inner structure using earthquakes: these 
generate seismic waves that travel through rock and other interior structures of our planet, 
which can then be sounded. Today, the Earth's interior has been completely mapped down 
to scales of a few hundred kilometres. 

The strongest earthquakes can even cause normal mode oscillations (e.g., see Montag- 
ner & Roult 2008), that is, the whole planet vibrates with its natural frequencies (also 
called eig en frequencies in theoretical analyses). These normal modes are most valuable 
in determining the deep interior structure of the Earth. Asteroseismology does just the 
same: it uses the frequencies of the normal modes of pulsating stars (that may be seen 
as "starquakes" ) as seismic waves. The eigenfrequencies of stellar models, that are de- 
pendent on their physical parameters and interior structures, are then matched to the 
observed ones. If a model succeeds in reproducing them, it is called a seismic model. 

The pulsation modes are waves in the stellar interior, just like the waves that musical 
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instruments resonate in. Popular articles often compare the frequency pattern of stellar 
oscillations to the sounds of musical instruments. The superposition of the normal mode 
frequencies in which a star oscillates can therefore be seen as its sound, although it 
generally does not involve frequencies the human ear is susceptible to. 

There is a large variety of normal modes that stars can pulsate in. The simplest are 
radial modes: the star periodically expands and shrinks, and its spherical symmetry is pre- 
served. The mathematical description of the displacement due to the oscillations results 
in differential equations of the Sturm-Liouville type that yields discrete eigensolutions: 
the radial mode frequencies of the given model. 

Pulsation in nonradial modes causes deviations from spherical symmetry: the star 
changes its shape. Mathematically, this no longer results in an eigenvalue problem of the 
Sturm-Liouville type, and a large number of possible oscillation modes originates. The 
eigenfunctions are proportional to spherical harmonics: 

y™(0, 0) = jV™pl m| (cos0)e^, (1) 

where 9 is the angle from the polar axis (colatitude), <fi is the longitude, P ; ' m ' is the 
associated Legendre polynomial, NJ 71 is a normalization constant and / and m are the 
spherical degree and azimuthal order of the oscillation. With 

pr(x) = (-ir(i-x*r /2 ^Pi(x) (2) 

we obtain 

P o °(cos0) = 1 
P°( C OS0) =COS0 Pi (cos 0) 

p 2 ° (cos 0) = ^ (3 cos 2 9-1) P\ (cos 9) = -3 cos 9 sin 9 

etc. etc. 

What does this mean in practice? Nonradial pulsation modes generate distortions on 
the stellar surface described by these spherical harmonics. The oscillations separate the 
stellar surface into expanding and receding as well as heating or cooling areas. A graphical 
example of these is shown in Fig. 2. 

Between the expanding and receding surface areas, no motion takes place. The lines 
along which this is the case are called the node lines; the nulls of the previously spec- 
ified associated Legendre polynomials specify the locations of these node lines in polar 
direction. The total number of node lines on the stellar surface is the spherical degree I, 
which must in all cases be larger than or equal to zero. The number of node lines that 
are intersected when travelling around the stellar equator is the azimuthal order m. I and 
m are the quantities that appear in the expressions for the eigenfunctions and for the 
associated Legendre polynomial given above, and they are used for the classification of 
the pulsation modes. 

Pulsation modes with m/0 are travelling waves (as can also be seen in the defining 
equation for spherical harmonics), and as they can run either with or against the rotation 



(3) 

= - sin 9 (4) 
P|(cos0) = 3sin 2 (5) 
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1=2, m=0 1=2, |=1 1=2, ^11 1=2 



Figure 2: Schematic description of the surface distortions produced by pulsation modes 
with spherical degrees < I < 2, arranged in the same way as the previous expressions 
for the associated Legendre polynomials. Whilst the outward moving areas of the star 
are coloured in dark grey, the light grey areas move inward, and vice versa. The pole and 
equator of the star are indicated. Adapted from Telting & Schrijvers (1997). 

of the star, m can lie in the interval [—/,£]. Modes with m = are called axisymmetric, 
modes with \m\ = I are named sectoral, and all other modes are referred to as tesseral 
modes. 

The third quantity needed to describe pulsation modes is the radial overtone k (some- 
times also denoted n in the literature), which is the number of nodes in the stellar interior. 
A mode that has no node in the interior is called a fundamental mode. A mode with one 
interior node is called the first overtone, modes with two interior nodes are the second 
overtone, etc. An accurate account of mode classification from the theoretical point of 
view is given by Deubner & Gough (1984). 

Historically, observationally, and inconsistently with the theoretical definition, radial 
overtone modes have also been called the first and second harmonics, respectively, and 
have been abbreviated with F for the fundamental, as well as 1H (or 10), 2H (or 20) 
etc. for the overtones. Radial pulsations can be seen as modes with £ = (remember 
that Pq(cos#) = 1); all other modes are nonradial oscillations. Modes with / = 1 are also 
called dipole modes; I = 2 modes are quadrupole modes. 

There are two major restoring forces for stellar oscillations that attempt to bring the 
star back in its equilibrium configuration: pressure and buoyancy (gravity). For radial 
motion, the gravitational force in a star increases during compression, so it would actually 
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accelerate, and not restore, the oscillation. Therefore, pressure must be the restoring force. 
On the other hand, for a predominantly transverse motion, gravity restores the motion 
through buoyancy, similar to what can be observed when throwing a stone into a pond. 
Therefore, aside from their identification with the pulsational quantum numbers k, I and 
to, nonradial pulsation modes are also classified into p pressure (p) modes and gravity 
(g) modes. These two sets of modes thus differ by the main direction of their motion 
(radial/transverse), and their frequencies. Pulsation modes with periods longer than that 
of the radial fundamental mode are usually g modes, whereas p modes have periods equal 
or shorter than that; radial pulsations are always p modes. 

The different modes are often labelled with their radial overtone number, e.g. a p 3 
mode is a pressure mode with three radial nodes, and a gg mode is a gravity mode with 
eight radial nodes. Modes with no interior nodes are fundamental modes, or f modes. 
Note that the f mode for I = 1 does not exist, as a dipole motion of the entire star would 
require a movement of the stellar centre of mass, which is physically impossible. 

The propagation of pulsation modes in the stellar interior is governed by two frequen- 
cies. One of these is the Lamb frequency Li, which is the inverse of the time needed to 
travel one horizontal wavelength at local sound speed. The other frequency describes at 
what rate a bubble of gas oscillates vertically around its equilibrium position at any given 
position inside a star; it is called the Brunt- Vaisala frequency N. These two quantities 
are defined as: 



where I is the spherical degree, c is the local velocity of sound, r is the radius, g 
is the local gravitational acceleration, p and p are local pressure and density in the 
unperturbed state, respectively, and 71 = (pdp/pdp)^ is the first adiabatic exponent. 

The Lamb and Brunt- Vaisala frequencies have the following implications: an oscillation 
with a frequency higher than both experiences pressure as the main restoring force in the 
corresponding part of the star. On the other hand, a vibration with a frequency lower 
than both L t and N is restored mostly by buoyancy. In other words, if we have a stellar 
oscillation with an angular frequency u>, it is a p mode wherever oj > Li,N, and it is a 
g mode wherever 00 < Li, N. In stellar interior regions where u lies between the Lamb 
and Brunt- Vaisala frequencies, the amplitude of the wave decreases exponentially with 
distance from the p and g mode propagation regions; such parts in the stellar interior are 
called evanescent regions. A propagation diagram aids the visualization of this discussion 



Whereas the Lamb frequency decreases monotonically towards the model's surface, 
the Brunt-Vaisala frequency shows a sharp peak near the model's centre and then rapidly 
drops to zero. This is because this model possesses a convective core, where N 2 < 0, and 
the spike in N is due to a region of chemical inhomogeneity. Over a range of stellar models, 
the behaviour of L\ in the interior is usually simple, whereas N may show considerable 
changes with evolutionary state and mass. However, N is always zero in the stellar centre. 

In Fig. 3, the g mode is confined to the innermost parts of the star. It is trapped in the 
interior, and therefore unlikely to be observed on the surface. The p mode is concentrated 
near the stellar surface and may be observable. The intermediate frequency mode shows 






(Fig. 3). 
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Figure 3: A propagation diagram for an evolved 2 M stellar model. The run of the Lamb 
(dashed line) and Brunt- Vaisala (full line) frequencies with respect to fractional stellar 
radius is shown. Some stellar pulsation modes are indicated with thin horizontal dashed 
lines, the circles are interior nodes. The lowest frequency oscillation (lowest dashed line) 
shown is a g$ mode; that with the highest frequency is a p$ mode. The oscillation with 
intermediate frequency is a mixed mode. Data kindly supplied by Patrick Lenz. 



remarkable behaviour: it has three nodes in the outer regions of the star, but also one 
node in the g mode propagation region. This particular mode is capable of tunnelling 
through its narrow evanescent region; it is a g mode in the deep interior, but a p mode 
closer to the surface. Such modes are called mixed modes. 

As mentioned before, stellar pulsation modes can be excited in certain parts of the 
stellar interior and they can propagate in some regions, whereas in other regions they are 
damped. The work integral is the energy gained by the pulsation mode averaged over 
one oscillation period. An evaluation of the work integral from the stellar centre to the 
surface is used to determine whether or not a given mode is globally excited in a stellar 
model. For excitation to occur, the exciting forces must overcome those of damping and 
the work integral W will be positive. 

The growth rate 7 = W/2nuI, where oj is the pulsation frequency and I is the mode 
inertia, parametrizes the increase of oscillation energy during a pulsation cycle, and also 
indicates how rapidly the amplitude of a given mode increases. The normalized growth 
rate (or stability parameter) 77, which is the ratio of the radius-integrated work available for 
excitation to the radius-integrated total work, is used to evaluate which pulsation modes 
are excited in the given model. If 77 > 0, a mode is driven and may reach observable 
amplitude; if rj = 1, a mode is driven in the entire stellar model, and if 77 = — 1, a mode 
is damped everywhere in the model. The most widely used application of this stability 
parameter is the comparison of the excited modal frequency ranges as determined from 
observations and those predicted by theory. 

Figure 4 shows how the frequency domains predicted to be excited change with the 
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evolution of a 1.8 M main sequence model. On the Zero-Age Main Sequence (ZAMS), 
only pure p modes of high radial overtone, and corresponding high frequencies, are excited 
in the model. Later on, some of these become stable again, but modes with lower overtone 
become excited. At the end of main sequence evolution, a large range of p modes, mixed 
modes and g modes is predicted to be excited, and the frequency spectrum becomes 
dense (and even denser as the model leaves the main sequence). The range of excited 
frequencies is an observable that, when compared with models, can give an estimate of 
the evolutionary state of the star. 
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Figure 4: Theoretically predicted / = 2 oscillation spectra of a 1.8 M & main sequence 
model, evolving from hotter to cooler effective temperature. Pulsation modes excited in 
this model are shown with filled circles, stable modes with open circles. The types of 
mode on the Zero-Age Main Sequence (ZAMS) are given. Note the g modes intruding 
into the p mode domain. 

The reason for the evolutionary change in the pulsation frequencies is that as the model 
evolves, its convective core, in which the g modes are trapped at the ZAMS, shrinks and 
the g mode frequencies increase. At the same time the envelope expands due to the 
increased energy generation in the contracting nuclear burning core, causing a decrease 
of the p mode frequencies. 

At some point, some p and g modes attain the same frequencies, and the modes begin 
to interact: they become mixed modes, with g mode character in the core and p mode 
character in the envelope (see also Fig. 3). The frequencies of these modes never reach 
exactly the same value; the modes just exchange physical properties. This effect is called 
avoided crossing or mode bumping (Aizenman, Smeyers & Weigert 1977). 

Of the individual modes, particular astrophysical potential is carried by the mode that 
originates as g\ on the ZAMS, as pointed out by Dziembowski & Pamyatnykh (1991). This 
mode is mostly trapped in the shrinking convective core (for stars sufficiently massive to 
possess a convective core), because the rapid change in mean molecular weight at its edge 
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causes a spike in the Brunt- Vaisala frequency (cf. Fig. 3). The frequency of this mode is 
thus dependent on the size of the convective core. 

The most important parameter determining the evolution of main sequence stars is 
the amount of nuclear fuel available. Stellar core convection can mix material from the 
radiative layer on top of the core into it, providing more nuclear fuel. This mixing is often 
named convective core overshooting and parametrized in theoretical models. Overshooting 
also decreases the gradient in the mean molecular weight at the edge of the core, which 
means that the frequency of the mode that is g\ on the ZAMS measures the convective core 
size, a most important quantity for astrophysics in general. To emphasize its sensitivity 
to the extent of stellar core convection, this mode has also been named g c mode. 

In the following, only p and g modes, and mixed modes of these types will be con- 
sidered. Other types of mode, such as r modes (torsional oscillations that may occur in 
rotating stars), g~ modes (convectively excited g modes in rotating stars), strange modes 
(showing up in calculations of highly nonadiabatic environments) or gravitational-wave 
w modes, will not be discussed as they have been of little practical importance for the 
seismic sounding of stars at the time of this writing. 

The frequencies of the p and g modes of pulsating stars depend strongly on their 
structure. However, when high radial overtones are considered, some simple relations 
between mode frequencies emerge. These are derived from asymptotic theory. The clas- 
sical reference on the subject is Tassoul (1980), the instructive is Gough (1986). In the 
high-overtone limit one finds 



where Au is the inverse sound travel time through the centre of the star, D is a 
frequency separation dependent on the stellar evolutionary state, Pq is the asymptotic 
period that is proportional to the integral of the Brunt- Vaisala frequency throughout the 
star, I and n are the spherical degree and radial overtone, respectively, and e and 5 are 
stellar structure parameters, respectively. 

These relations have important consequences: low-degree p modes of consecutive high 
radial overtones of the same spherical degree are equally spaced in frequency, whereas 
low-degree g modes of consecutive high radial overtones of the same spherical degree are 
equally spaced in period. Furthermore, if the parameter Dq were zero, Eq. 7 indicates 
that the frequencies of high-overtone p modes of even and odd degrees would be the same, 
respectively, and odd-degree modes would have frequencies intermediate between those 
of even degree. In realistic stellar models, these relations hold approximately, but not 
exactly. 

In a nonrotating star, the frequencies of modes with m^O are the same as those of 
the m = modes. However, as the m^O modes are travelling waves, their frequencies 
separate in the observer's frame when looking at a rotating star: the mode moving with 
rotation appears at higher frequency, the mode moving against rotation appears at lower 





for p modes, whereas for g modes 




(8) 
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frequency and the frequency difference to their nonrotating value is m times the rotation 
frequency (e.g., see Cox 1984). This effect is called rotational frequency splitting, and is 
one basic tool of asteroseismology: if such splittings are observed, the rotation frequency 
of the star can be determined. 

Unfortunately, reality is not quite as simple as that. The Coriolis force acts on the 
travelling waves and modifies their frequencies. In addition, the m^O modes cause tidal 
bulges, on which centrifugal forces act. Therefore, the frequencies of these modes are 
often expressed as 

Uk,i, m = ^k,ifl + rn(l-C K i)VL + m 2 ^— (9) 

CO 

in case of (moderately) slow stellar rotation (Dziembowski & Goode 1992). iOkj >m is the 
observed frequency of the mode k, I, m, Q is the stellar rotation frequency, and Ck,i and Dk,i 
are constants that describe the effects of the Coriolis and centrifugal forces, respectively; 
they are also called Ledoux constants. These constants are usually determined from stellar 
model calculations. The rotational splitting constant also approaches asymptotic 
values for high radial overtones. For g modes, C^,i becomes (1(1 + l)) -1 , and for p modes 
Ck,i = in the asymptotic limit. 

It should be made clear that the previous formula is only an approximation for the case 
that rotation can be treated as a perturbation to the equilibrium state: stellar rotation 
distorts the spherical shape of a star. As a consequence, the individual modes can no 
longer be described with single spherical harmonics. For instance, radial modes receive 
contamination of I — 2 and other modes with even /. Vice versa, I = 2 modes obtain some 
1 — 0,1 — 4, etc. contributions. In the previous formula this means that for rapid rotation 
higher order terms may be added, but to arrive at reliable results, two-dimensional numer- 
ical calculations are really required. Rotational mode coupling can affect the properties 
of oscillation modes with close frequencies (e.g., see (Daszyhska-Daszkiewicz et al. 2002). 
Furthermore, the rotational distortion causes the stellar temperature to increase at the 
flattened poles and to decrease at the equatorial bulge. It follows that asteroseismology 
of rapidly rotating stars has an additional degree of complexity. 

Even more complexity to the theoretical treatment of stellar pulsation is added by 
the presence of a magnetic field. A weak magnetic field would generate a second-order 
perturbation to the pulsation frequencies, just like the centrifugal force in case of rotation, 
but with opposite sign, i.e. the observed oscillation frequencies would increase with respect 
to the non- magnetic case (e.g., Jones et al. 1989). 

In the presence of a strong magnetic field, the effect would be more severe: if the 
field was oblique to the rotation axis (just like the Earth's magnetic field!), the pulsation 
axis would align with the magnetic axis, and no longer with the rotation axis, as implic- 
itly assumed so far. This means that an observer would see each pulsation mode at a 
varying angle over the stellar rotation period, creating amplitude and phase variations 
with exactly that period, and therefore separable from the effects of rotational m-mode 
splitting. A star oscillating in this way is called oblique pulsator (Kurtz 1982). We note 
for completeness that the most general case of rotational splitting is the nonaligned pul- 
sator (Pesnell 1985), where (21 + l) 2 components of a given mode appear, with frequency 
separations proportional to the rotation period and to the angle between the rotation and 
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pulsation axes. 

To summarize, different stellar pulsation modes propagate in different interior regions, 
and their energy within those regions is not equally distributed. Each single pulsation 
mode has a different cavity, and its oscillation frequency is determined by the physical 
conditions in its cavity. This means that different modes are sensitive to the physical 
conditions in different parts of the stellar interior. Some modes teach us more about 
stellar envelopes, whereas other modes tell us about the deep interior. The more modes 
of different type are detected in a given star, the more complete our knowledge about its 
inner structure can become. 

Some stars do us the favour to oscillate in many of these radial and/or nonradial 
modes simultaneously. Interior structure models of the stars can then be refined by 
measuring the oscillation frequencies of these stars, identifying them with their pulsational 
quantum numbers, and by reproducing all these with stellar models. This method is 
very sensitive because stellar oscillation frequencies can be measured to extremely high 
precision, up to one part in 10 9 (Kepler et al. 2005). In the following it will be described 
what observables are available to base asteroseismic models on, how the measurements 
can be interpreted, and how observations and theory are used to develop methods for 
asteroseismic interpretations. 

From the telescope to a seismic model 
Basic methods for analysing asteroseismic data 

Because stellar oscillations generate motions and temperature variations on the surface, 
they result in observable variability. The interplay of these variations causes light, radial 
velocity and line profile changes. Pulsating stars can thus be studied both photometrically 
and spectroscopically, via time series measurements. 

These time series are subjected to frequency analysis, meaning that the constituent 
signals are extracted from the data. In many cases, this is done by harmonic analysis, 
transforming the time series into frequency/amplitude space, e.g. by using the Discrete 
Fourier Transformation 



of the input function x(tk), corresponding to the time series of the measurements. 
Periodograms, amplitude or power spectra are means of visualizing the results; an example 
is given in Fig. 5. 

The amplitude spectrum in Fig. 5 can be used to estimate the frequencies of the 
dominant signals in the time series. In many cases, the analysis is carried forward by 
fitting sinusoids to the data, determining and optimizing their frequencies, amplitudes 
and phases, often by least squares methods. It has also become common practice to 
subtract this optimized fit from the data and to compute periodograms of the residuals, a 
procedure called prewhitening. This process is repeated until no more significant signals 
can be found in the data (a most delicate decision!). Depending on the specific case and 
requirements of the data set, a large number of alternative frequency analysis methods can 
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Figure 5: The light variations (upper two panels) of a pulsating star and the corresponding 
Fourier amplitude spectrum. The complicated beating in the light curve is reflected by 
the presence of several signals in the periodogram; the strongest are labelled with arrows. 



also be applied, such as phase-dispersion minimization, autocorrelation methods, wavelet 
analysis etc. Care must be taken to keep the limitations of the data sets and of the applied 
methods in mind. Periods no longer than the data set itself can be reliably determined, 
and adjacent frequencies spaced no more than the inverse length of the data set can be 
resolved. Nonsinusoidal signals cause harmonics and combination frequencies in Fourier 
amplitude spectra that must not be mistaken for independent mode frequencies. No 
frequencies shorter than the inverse of twice the median distance between consecutive 
data points can be unambiguously detected; the highest frequency that can be retrieved 
in a given set of measurements is also called its Nyquist frequency. 

Having determined the frequencies characterizing the stellar variability, the next step 
is their interpretation. Measurements of distant stars have an important limitation: non- 
radial oscillations create patterns of brighter and fainter, approaching and receding areas 
on the stellar surface. However, as a distant observer can usually not resolve the stellar 
surface, she or he can only measure the joint effect of the pulsations in light and radial 
velocity. As a consequence, the effects of oscillations with high spherical degree average 
out in disk-integrated measurements, and their observed amplitudes are reduced with 
respect to the intrinsic value. This effect is called geometric cancellation (Dziembowski 
1977). Calculations show that the amplitude drops as ~ \j\fl for high /. In ground based 
observational studies it is mostly assumed (and confirmed) that only modes with / < 2 are 
observed in light and radial velocity, with a few exceptions of I up to 4. Radial velocity 
measurements are somewhat more sensitive to higher / than photometric observations. 
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Once the oscillation frequencies of a given star have been determined, how can they be 
made asteroseismic use of? As an example, high-precision radial velocity measurements 
of the close-by star a Centauri A showed the presence of solar-like oscillations; the power 
spectrum (square of amplitude vs. frequency) of these data is shown in Fig. 6. 



0.2 




2.4 2.6 
Frequency (mHz) 



Figure 6: High-overtone p mode spectrum of the solar-like oscillator a Centauri A (Bed- 
ding et al. 2004). The vertical dotted lines are separated by 106.2 /iHz. 



This graph contains a series of maxima equally spaced in frequency: a high-overtone 
p mode spectrum, as predicted by asymptotic theory. The mean frequency spacing is 
106.2 /iHz. However, it is obvious that only every other of the strongest peaks conforms to 
this spacing; there are other signals in between. The signals halfway between the vertical 
lines in Fig. 6, denoting modes of spherical degree /, are a consequence of Eq. 7: these 
are modes with I ± 1. Given the effects of geometrical cancellation, it is straightforward 
to assume that these are the modes of lowest /, viz. / = and 1. 

However, there is more information present. A close look at Fig. 6 reveals that many 
of the strongest peaks seem to be split. Again, the explanation for this finding lies within 
Eq. 7: the close neighbours are modes with 1 + 2. These frequencies no longer have the 
same values as those with degree / because the stellar interior has structure. 

As stars evolve on the main sequence, their nuclear burning cores shrink, increase 
in density, and change in chemical composition as evolution progresses. This alters the 
acoustic sound speed in the core and is reflected in the frequency differences between 
modes of degrees /, 1+2, also called the small frequency separation 5u m . On the other hand, 
evolution causes expansion of the outer regions of stars that become more tenuous, which 
increases the sound travel time through the star. This is measurable via the frequency 
difference between consecutive radial overtones and is called the large frequency separation 
Au . 

The large and the small separations can be computed for a range of theoretical stellar 
models. It turns out that a plot of 5u vs. Az/ allows an unambiguous determination of 
stellar mass and evolutionary state, and this method works particularly well for models 
with parameters similar to our Sun (main sequence, M £ 1.3M Q ). This diagnostic is 
called an asteroseismic HR Diagram (Christensen-Dalsgaard 1988). 

Another important tool to analyse high-overtone p mode pulsation spectra is the 
Echelle Diagram. This diagram plots the oscillation frequencies versus their modulus 
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with respect to the large separation. Figure 7 shows the Echelle Diagram for a Cen A 
constructed from the frequencies of the signals apparent in Fig. 6. 
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Figure 7: Echelle diagram of the frequencies determined from Fig. 6 (values from Bedding 
et al. 2004). Full dots represent radial modes, triangles mark dipole modes, plus signs 
stand for quadrupole modes and diamonds are for / = 3 modes. 



The frequencies fall onto four distinct ridges, corresponding to modes of the same 
spherical degree /. In this way, I can be identified. From left to right, the ridges correspond 
to I = 2, 0, 3 and 1, respectively, and it can be seen that, within the errors, 8vi% = 5/3<5z/o2, 
consistent with Eq. 7. The scatter in this diagram is due to a combination of the temporal 
resolution of the data and, more importantly, the finite lifetimes of the stochastically 
excited modes. 

The example of a Cen A shows that once a sufficient number of intrinsic pulsation 
frequencies of a given star has been determined, one may identify their spherical degrees by 
pattern recognition. This method is also applicable to the oscillation spectra of pulsating 
white dwarf stars, as exemplified in Fig. 8. 

The series of peaks discernible in Fig. 8 does not form a pattern of equally spaced 
frequencies, but of equally spaced periods. This is the signature of a high overtone g 
mode pulsation spectrum (Eq. 8). In this case, all the strong peaks are roughly aligned 
with the vertical lines; there are no strong modes in between. Some of the expected signals 
are missing, but in most cases only apparently: they are just much weaker in amplitude. 
A mean period spacing of 39.5 seconds has been determined from this analysis, and 
the obvious identification of the strongest modes is I — 1. Equation 8 then results in 
Po = 55.9 s. For pulsating white dwarf stars, Pq is a measure of the stellar mass, which 
was consequently determined. 

Furthermore, several of the strongest peaks in Fig. 8 are split into triplets. These 
triplets are equally spaced in frequency, and they are the signature of rotationally split m 
modes (Eq. 9). The rotation period of GD 358 could therefore also be determined, just 
from an inspection of the power spectrum of its photometric time series, and basic appli- 
cation of theory. The fact that the m-mode splitting only results in triplets strengthens 
the previous identification of the spherical degree as / = 1. The relative amplitudes within 
those triplets are dependent on the inclination of the stellar pulsation axis (Eq. 4), which 
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Figure 8: High-overtone g mode spectrum of the pulsating DB white dwarf star GD 358 
(Winget et al. 1994). The vertical dotted lines denote periods with an equal spacing of 
39.5 seconds. 



may then be determined, but in reality there are other, presently unknown, effects that 
modify the relative multiplet amplitudes. 

There is even more information present. It is noticeable that the peaks in the amplitude 
spectrum do not perfectly conform to their asymptotically predicted locations. This is a 
sign of mode trapping. White dwarf stars consist of a degenerate core, with subsequent 
outer layers of different chemical elements. The transition regions between these layers 
create spikes in the Brunt- Vaisala frequency. The pulsation modes prefer to place their 
radial nodes in these transition regions, and be standing waves to both sides of the nodes. 
This modifies their frequencies compared to the case of homogeneous interior structure 
and gives rise to the observed deviations from equal period spacing. 

The observant reader will have noticed three apparent inconsistencies in the previous 
paragraphs. How can two stars with similar oscillation frequencies be high-overtone p 
and g mode pulsators, respectively? How does one know that 106.2 /zHz is the mean 
frequency spacing for a Cen A, and not half this value? Why can the claimed small 
frequency separation not be the effect of rotation? 

This is because some additional constraints are available that help in the interpretation 
of pulsational mode spectra. As the p mode pulsation periods depend on the sound travel 
time through the star, they must be related to its size, or more precisely, to its mean 
density. The pulsation constant 

Q = pM, (ii) 
V p 

where P is the pulsation period and p is the mean stellar density, is a useful indicator 
of what type of mode one sees in a given star. Over the whole HR diagram, the Q value for 
the radial fundamental mode is between 0.03 — 0.04d; for the Sun, it is 0.033d (0.8 hr). 
As mentioned before, radial modes can only be p modes, (pure) p modes always have 
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frequencies of the same order or larger than that of the radial fundamental, and (pure) g 
modes always have frequencies lower than that. 

Because GD 358 is a white dwarf star and therefore has high mean density, its radial 
fundamental mode period would be of the order of 4 s. On the other hand, a Cen A, 
a little more massive and more evolved than our Sun, has a radial fundamental mode 
period of about 1 hr. Therefore, the similar 5-10 minute pulsation periods of the two 
stars correspond to completely different types of mode. 

The asymptotic p mode frequency separation also relates to the stellar mean density 



where Az/ 0iQ = 135yuHz. Given the knowledge on the mass and evolutionary state of 
a Cen A, it immediately be inferred that 106.2 /iHz must be the large frequency spacing. 
Finally, the rotation period of a Cen A is much too long to generate m-mode splitting 
with the observed <5z/ 2. 

Apart from these examples, there is another clue towards the nature of observed pul- 
sation modes just from observed frequency spectra. Radial modes are global oscillations 
that may reveal themselves because the period ratios of consecutive overtones are well 
known. On the main sequence, the period ratio of the radial fundamental mode and the 
first overtone is around 0.773, and the period ratio of the first to second overtone radial 
modes is 0.810. For more evolved stars, such as the 5 Cephei stars, these ratios change to 
0.705 and 0.805, respectively. Such period ratios are therefore suggestive to correspond 
to radial modes. 

Methods for mode identification 

It is now clear that besides the pulsational mode spectra of stars themselves, the incorpo- 
ration of other constraints is useful for their interpretation. This is particularly important 
for oscillation spectra that do not show obvious imprints of the underlying modes, like 
those of stars pulsating in high radial overtones. Examples are pulsations of low radial 
overtone or stars rotating so rapidly that the rotational splitting is of the same order as the 
frequency spacing of consecutive radial overtones of the same I or of modes of the same k, 
but different I. In practice, most observed stellar oscillation spectra are incomplete, either 
because the star chooses so or because the observations do not have sufficient sensitivity, 
which makes it even more difficult to recognize patterns in the observed frequencies or 
periods and to type the modes accordingly. Therefore, mode identification methods have 
been developed, of which a variety is available. 

The first method uses photometric data only. The flux change for nonradial pulsation 
in the linear regime can be expressed as 



Am(X,t) = -1.086ePrW((7i + T 2 ) cos(wt + Vt) + (T 3 + T 4 + T 5 )cos(cut)) (13) 



(Watson 1988), where Am(A, t) is the time and wavelength dependent magnitude vari- 
ation of an oscillation, — 1.086e is an amplitude parameter transformed from fluxes to 




(12) 
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magnitudes, P ; ' m ' is the associated Legendre polynomial, fiQ is the cosine of the inclina- 
tion of the stellar pulsation axis with respect to the observer, oj is the angular pulsation 
frequency, t is time and i{jt is the phase lag between the changes in temperature and local 
geometry (mostly originating in convection zones near the stellar surface). The term T\ 
is the local temperature change on the surface, T 2 is the temperature-dependent limb 
darkening variation, T 3 is the local geometry change on the stellar surface, T 4 is the local 
surface pressure change and T 5 is the gravity-dependent limb darkening variation. 

The Tj terms can be determined for different types of pulsator from theoretical model 
atmospheres, and the observables best suited to reveal the types of mode present can 
also be deduced. These would for instance be the photometric amplitude ratios or phase 
shifts between different filter passbands, but also the optimal passbands themselves. An 
example of such a mode identification is shown in Fig. 9. 
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Figure 9: Identification of the four strongest pulsation modes of the (3 Cephei star 12 Lac- 
ertae from multicolour photometry (taken from Handler et al. 2006). The amplitudes are 
normalized to unity in the ultraviolet and compared with theoretical predictions. The full 
lines are for I = 0, the dashed lines for / = 1, the dashed-dotted lines for / = 2, the dotted 
lines for / = 3, and the dashed-triple dotted lines for I = 4. The modes investigated in 
the two left panels are 1 = 1, the next is / = 2 and the rightmost one is / = 0. 

From photometry alone, only the spherical degree of a given pulsation mode can be 
identified. The method can be supported by adding radial velocity measurements, which 
increases its sensitivity (Daszyhska-Daszkiewicz, Dziembowski, & Pamyatnykh 2005), but 
still does not supply a determination of the azimuthal order. To this end, high resolution 
spectroscopy must be invoked. 

The lines in a stellar spectrum are broadened by rotation through the Doppler effect: 
the intrinsic line profile is blueshifted on the parts of the stellar surface approaching the 
observer, and redshifted on the areas moving away. The effect is strongest on the stellar 
limb and decreases towards the centre. As a consequence, a rotationally broadened line 
profile contains spatial information of the stellar surface. For instance, in a dark starspot 
flux is missing, and a spike will move through the line profile as the spot rotates over 
the visible disk. The method to reconstruct stellar surface structure from line profile 
variations is called Doppler Imaging (Vogt & Penrod 1983). 

Apart from radial velocity changes, stellar oscillations also cause variations in rota- 
tionally broadened line profiles. The areas on the stellar surface that have an additional 
approaching component to their motion with respect to the observer have their contri- 
bution to the line profile blueshifted, whereas the receding parts are redshifted by the 
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corresponding amount of pulsational Doppler velocity. The net result of all these mo- 
tions is bumps travelling through the line profile, and their shapes are governed by the 
oscillation mode causing them (e.g., see Telting 2003 for a review). Examining stellar 
line profiles, pulsation modes of up to I ~ 20 can be observed and identified, a vast ex- 
tension in / compared to photometric and radial velocity techniques. Some examples of 
pulsational line profile variations are shown in Fig. 10. 




Figure 10: Line profile variations due to stellar pulsation. The upper parts of the graph 
show the shape of the oscillation mode on the surface, whereas the thin lines in the 
lower halves represent the unperturbed rotationally broadened line profile, and the thick 
lines are the superpositions with the pulsation. Each individual mode (from left to right: 
I — 4, m — 0; I — 5, \m\ = 3; I = \m\ = 7) generates a different distortion of the line profile. 
Adapted from Telting & Schrijvers (1997). 



The task now is to extract the correct values of I and m from the observed variations. 
The principle is to fit the theoretically calculated 3-D velocity field to the observed line 
profiles, and a wide range of spectroscopic mode identification methods is available. Some 
of the most commonly used are the Moment method (most suitable for low I), the pixel- 
by-pixel method, the Fourier Parameter Fit method, or Doppler reconstruction of the 
surface distortions. High-resolution spectroscopy is better suited for the determination 
of m rather than / which makes it complementary to photometric and radial velocity 
methods. 

When resorting to photometric or spectroscopic methods, it is not required to arrive 
at unique identifications of all observed pulsation modes in each given star. What is 
needed is the secure identification of a sufficient number of modes to rule out all possible 
alternative interpretations. An example will be presented later. 
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Asteroseismic modelling 

The observed pulsation frequencies and their identifications are then matched to theoreti- 
cal models (see Kawaler, this volume, for details). These would ideally be full evolutionary- 
models with pulsation codes operating on them, although a few codes based on envelope 
models are still in use. 

Most pulsation codes use the linear approximation: the oscillations are treated as linear 
perturbations around the equilibrium state. This allows the evaluation of the excitation of 
oscillation modes and thus the computation of theoretical domains of pulsation in the HR 
diagram. Nonlinear computations, that would allow predictions of oscillation amplitudes, 
are still rather the exception than the rule because they are, even today, expensive in 
terms of computing time, numerical hydrodynamical simulations. 

Many stellar pulsation codes employ the adiabatic approximation (sometimes called 
isentropic approximation) to compute oscillation frequencies, which is the assumption 
that no energy exchange between the oscillating element and its surroundings takes place. 
Other codes perform nonadiabatic frequency computations. A wide variety of stellar 
oscillation codes is available, and most theory groups use their own routines, optimized 
for the application to their objects of main interest. 

There are several strategies to find seismic models from observations. Some compare 
the observed oscillation frequencies with those of a grid of stellar models and perform 
automatic matching between them. This is computationally extensive, which means that 
supercomputers or parallel processing are invoked, or that intelligent optimization meth- 
ods such as genetic algorithms are employed, or both. 

Other strategies start with first imposing observational constraints. Besides the ob- 
served oscillation frequencies themselves and the identification of the underlying pulsation 
modes, these would often be estimates of the objects' positions in the HR diagram. As 
an example, these are shown for the (3 Cephei star v Eridani in Fig. 11. 

The detection of a radial mode in the frequency spectrum is an asset for asteroseismic 
studies of this star: there are only three possibilities for the value of its mean density. 
These depend on whether its frequency corresponds to the fundamental, the first or the 
second radial overtone (Eq. 11). A comparison of the observed position of the star in 
the HR Diagram (right panel of Fig. 11) with its error bars leads to the rejection of the 
second overtone hypothesis, and to the elimination of models with masses below 8.5 M . 

Now the / = 1 modes come into play: moving along the lines of constant mean density 
in the HR diagram, a comparison between their observed and theoretically predicted 
frequencies can be made. This is done on the left side in Fig. 12. Given the uncertainties 
and assumptions in the asteroseismic model construction, all / = 1 modes are reproduced 
by models between 9.1 — 9.8 M that have the radial mode as the fundamental. However, 
the frequency of the highest-overtone dipole mode cannot be explained within the errors 
for models that assume the first overtone for the radial mode. This interpretation can 
therefore also be rejected. The modal assignment to the observed frequencies is now 
unambiguous, and the range of models to be explored for seismic fitting is severely reduced. 

The reason why this way of mode identification was successful is that the g\ mode and 
its respective p mode neighbour are in the process of avoided crossing in the parameter 
space of interest. Small changes in the evolutionary states of the models depending on 
mass therefore lead to significant changes in their frequencies. The observed / = 2 mode 
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Figure 11: Left: schematic oscillation spectrum of v Eridani. The numbers on top of each 
mode (group) are their I identifications, consistent in photometry and spectroscopy (De 
Ridder et al. 2004). Right: a plot of the star's position in the theoretical HR Diagram (star 
symbol) with its error bars and lines of equal mean density for the observed radial mode 
periods assuming they the fundamental, first and second overtones, respectively (thick 
lines). Some model evolutionary tracks labelled with their masses, the Zero- Age Main 
Sequence, and the borders of the /3 Cephei (dashed-dotted line) and SPB star (dashed 
lines) instability strips are also shown. 



has been excluded from the mode identification procedure as all models in the parameter 
domain under consideration reproduce its frequency correctly: it is nearly a pure p mode 
and its frequency hardly changes in models of the same mean density. 

In this restricted parameter space, an interesting observation can be made: the ob- 
served / = 1 m-mode splittings do not agree with those predicted by uniformly rotating 
models (Fig. 12, right-hand side). Fitting the rotational splitting for the g\ mode, that 
samples the deep stellar interior, results in predicted splitting about 30% larger than that 
observed for the pi mode, which is concentrated closer to the surface. This means that 
the star's rotation rate increases towards its interior, as predicted by theory (Talon et al. 
1997). In addition, as the I — 1, g\ mode frequency is sensitive to the size of the convective 
core (it is the g c mode!), a constraint on the convective core overshooting parameter was 
also obtained (Pamyatnykh, Handler, & Dziembowski 2004) 

This example shows the potential of asteroseismology, even if only few pulsation modes 
are available. Radial pulsations are most valuable in the identification process as they 
immediately provide accurate constraints on the stellar mean density that can become 
unambiguous if the radial overtone of the mode can be inferred, supported by other 
estimates such as on effective temperature and luminosity. Theory and observations work 
hand in hand. Once a unique identification of the normal modes is achieved, firm results 
on stellar structure can be obtained, even given uncertainties in the modelling procedure. 
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Figure 12: Left: a figure showing the match of the I — 1, m — modes of v Eri for models 
of different mass but same mean density. The full horizontal lines are the observed 
frequencies, the dotted lines are theoretical model frequencies. Upper panel: assuming 
that the radial mode is the fundamental. The vertical dashed-dotted lines show the 
mass range in which the first g mode and the first p mode fit the observed frequencies 
(9.1 — 9.8 M Q ). Lower panel: the same, but assuming the radial mode is the first overtone. 
Here it would be the first g mode and the second p mode that give an acceptable fit 
between 10.0 — 10.8 M . However, the observed I = 1 mode with highest frequency is not 
compatible with the more massive models. Right: a comparison of the rotational splittings 
of a rigidly rotating model, with a rotation rate chosen to fit the observed / = l,gi triplet, 
with the observations. The observed splitting of the I = l,pi triplet is not reproduced 
with this assumption, demonstrating the presence of differential interior rotation. 



In the case of v Eridani, the most severe problem is that the / = l,p2 mode is hard to 
be excited and accurately matched in frequency given the observed effective temperature 
of the star and its chemical surface composition. Only changes to the input physics would 
solve this problem. This is just the goal of asteroseismology: to improve our knowledge 
about stellar physics. Such an improvement can only be achieved using observational 
results that present models cannot account for! 



Applications 

Pulsating white and pre-white dwarf stars 

The first pulsators that could be studied asteroseismically were white dwarf and pre-white 
dwarf stars, pulsating in high-overtone low-degree g modes. A remarkable initial result 
was the theoretical prediction of a then new class of pulsating white dwarfs of spectral 
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type DB, and their subsequent observational confirmation (Winget et al. 1982). Upon 
the realization that some pulsating white dwarf stars have very complicated frequency 
spectra, and having a well-developed theory of white dwarf pulsation available, the main 
obstacle to extracting correct mode frequencies from observational data was the daily 
interruptions of the measurements by sunlight. Figure 13 explains why. 




J I I I I I I I I I L 



0.5 1 

Time (days) 

Figure 13: Simulation of observations from a single astronomical site. The filled circles 
represent the measurements, including noise. The full line is a fit with the correct fre- 
quency present in the data, and the dotted line is a fit with a frequency different by 1 
cycle per day. The two fits represent the data almost equally well. However, they would 
be completely out of phase, and therefore easily separable, if measurements were available 
in between the present data, from a site 180 degrees different in geographical longitude. 
This led to the setup of multisite telescope networks. 

The interruptions cause ambiguities in the determination of the frequencies of the sig- 
nals present in the data: a frequency different by one cycle per sidereal day from the real 
oscillation frequency would generate a fit of comparable quality. In the presence of com- 
plicated variability and most notably for signals of low signal-to-noise this could lead to 
erroneous frequency determinations. Any seismic model based on incorrect observational 
input is misleading. 

The solution to this problem is to avoid, or at least minimize, daytime gaps in time 
resolved measurements. This can be accomplished by concerted observational efforts, 
involving interested colleagues over the whole globe, passing on the asteroseismic target 
from one observatory to the next. The best known of these collaborations is the Whole 
Earth Telescope (WET, Nather et al. 1990), invented for the study of pulsating white 
dwarf stars. 

One of the first WET runs was devoted to the prototypical pulsating pre-white dwarf 
star PG 1159-035 = GW Vir. It resulted in the detection of over 100 g-mode frequencies of 
I = 1 and 2 modes of high radial overtone, leading to precise determinations of the stellar 
mass and rotation period, an asteroseismic detection of compositional stratification, and 
an upper limit to the magnetic field strength (Winget et al. 1991). Subsequent WET 
observations of the prototypical pulsating DB white dwarf star GD 358 = V777 Her 
showed a mode spectrum dominated by high-overtone / = 1 g modes (cf. Fig. 8), resulting 
in determinations of its total and surface Helium layer mass, luminosity and rotation rate 
(Winget et al. 1994). 
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The evolution of white dwarf stars is dominated by cooling, at (nearly) constant radius. 
As they cool, they pass through a number of instability strips. Evidence is that all white 
dwarf stars located in such instability domains in the HR diagram do pulsate. This has 
an important consequence: the interior structures of the pulsators must be representative 
of all white dwarf stars. Thus asteroseismic results for white dwarf pulsators can be 
extended to the stars that do not oscillate. 

Cooling of pulsating white dwarf stars changes their oscillation periods, and the rate 
of period change is directly related to their energy loss. The hottest DB pulsators have 
reasonably high neutrino emission rates, and their evolutionary period changes may be 
able to tell us whether their neutrino emission is compatible with the standard model of 
particle physics (Winget et al. 2004). Measurements to detect such a period change are 
ongoing. On the other hand, the period changes of DA pulsators (some of which are the 
most stable optical clocks in the Universe, Kepler et al. 2005) could reveal the masses of 
axions, if the latter existed (Kim, Montgomery, & Winget 2006). 

Asteroseismology of white dwarf pulsators does not only allow to detect stratification 
in their chemical profiles near the surface, it also gives evidence of their core composition. 
This, in turn, is dependent on their history of evolution on the Asymptotic Giant Branch 
(AGB), and can be used to obtain constraints on the nuclear reaction rates in AGB stars. 
Present results (Metcalfe 2005) indicate consistency with evolutionary models. 

As white dwarf stars cool, their cores become crystallized. Being composed mainly of 
carbon and oxygen, such cores can be seen as giant diamonds! Massive white dwarf stars 
begin to crystallize when still in the DAV (ZZ Ceti) instability strip, and asteroseismic 
investigations of one such massive pulsator have proven substantial crystallization in its 
interior (Metcalfe, Montgomery & Kanaan 2004, Brassard & Fontaine 2005). 

The light curve shapes of pulsating white dwarf stars of spectral types DB and DA are 
often nonsinusoidal. The nonlinearities originate in their convection zones, that cannot 
instantly adjust to the pulsational variations (the g mode pulsations of pulsating white 
dwarf stars are almost exclusively due to temperature changes; Robinson, Kepler, & 
Nather 1982). As the light curve shapes of such pulsators depend on the thermal response 
time of the convection zone, the latter parameter can be determined from nonlinear light 
curve fits (Montgomery 2005). 

As a final example, there are pulsating white dwarf stars in mass accreting close binary 
systems. If the mass transfer rate is in a certain range, the surface temperature of the 
accreting white dwarf places it in an instability strip. A handful of such oscillators is 
known to date (Mukadam et al. 2007), but attempts at asteroseismology have proven 
difficult due to low amplitudes and unstable mode spectra. 

Among all classes of pulsating star (aside from the Sun itself), asteroseismology of 
pulsating white dwarf stars is certainly in the most advanced state. A new class of such 
oscillators has recently been proposed, hot DQ stars with carbon-dominated atmospheres 
and temperatures similar to that of the DB pulsators (Montgomery et al. 2007, Dufour 
et al. 2009). There is little doubt remaining that the variability of these stars is due to 
pulsation. Theory predicts yet another new type of white dwarf oscillator, DA white dwarf 
stars somewhat hotter than the DB pulsators, but observational searches for them have so 
far been inconclusive (Kurtz et al. 2008). We refer to Montgomery (2009), and references 
therein, for more information on asteroseismology of pulsating white dwarf stars. 
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Delta Scuti stars 

At about the same time when the necessity for worldwide observing efforts for pulsating 
white dwarf stars was realized, the same conclusion was reached for 5 Scuti stars, some 
of which also exhibit complex oscillation spectra composed of g, p and mixed modes of 
low overtone. The Delta Scuti Network, founded over 25 years ago, was the first multisite 
observing collaboration for these stars, followed by a number of others such as STEPHI 
or STACC; a few WET runs were also devoted to S Scuti stars. 

The asteroseismic potential of 5 Scuti stars is enormous, but could so far not be fully 
exploited. Part of the reason is visible in Fig. 4: the pure p mode spectrum on the 
ZAMS mainly allows a determination of the stellar mean density. When the scientifically 
more interesting mixed modes appear, the frequency spectrum is fairly dense and requires 
a large amount of data and long time base to be resolved observationally. However, 
the real stars seldom co-operate in showing many of the potentially excited modes at 
observable amplitude, inhibiting mode identification by pattern recognition. In addition, 
mode amplitudes are often small, counteracting reliable identifications of many modes by 
applying the methods discussed before. Even though dozens or even hundreds of pulsation 
modes have been detected in some 5 Scuti stars, little has been learnt on their interior 
structures from asteroseismology so far. 

Pulsational amplitude limitation of S Scuti stars is a major problem for theory: what 
makes the star excite only certain modes, which modes would these be and what de- 
termines their amplitudes? Observational evidence suggests that in evolved stars mode 
trapping is (part of) the answer, and that oscillations with frequencies around those of 
radial modes are preferentially excited (Breger, Lenz, & Pamyatnykh, 2009). 

For slowly rotating 5 Scuti stars gravitational settling and radiative levitation give 
rise to chemical surface peculiarities and are believed to deplete the pulsational driving 
regions. Consequently, Am and Ap stars are not expected to pulsate, although a few of 
them do (Kurtz 2000). This also means that many 5 Scuti stars rotate rapidly, which 
requires special calculations to extract information from their distorted pulsation modes, 
a field that has made considerable progress in the recent past (Reese et al. 2009). 

Some pre-main sequence stars cross the 5 Scuti instability strip on their way to the 
ZAMS, and consequently pulsate. The interior structures of these stars are fairly simple, 
which may make them more accessible to asteroseismic investigation compared to their 
main sequence counterparts. The oscillation spectra of pre-main sequence and main se- 
quence 5 Scuti stars in the same position in the HR Diagram are predicted to be different, 
which may allow an asteroseismic separation (Suran et al. 2001). 

Slowly Pulsating B and Gamma Doradus stars 

These two classes of high-overtone g mode pulsator, although well separated in effective 
temperature, share most of their asteroseismic characteristics. They also share the prob- 
lems with respect to observations, caused by their long periods: resolving their oscillation 
spectra requires measurements over a long time baseline, possibly many years. It is there- 
fore no surprise that most of the known members of these two groups of pulsator were 
discovered with the HIPPARCOS (High Precision PARallax Collecting Satellite) mission 
from its data set spanning over three years (Waelkens et al. 1998, Handler 1999). 
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As in pulsating white dwarf stars, effects of inhomogeneities in stellar structure would 
manifest themselves in mode trapping and thus in oscillatory behaviour in the g mode 
periods (Miglio et al. 2008). The dominant inhomogeneity is the change in mean molecular 
weight at the edge of the convective core, whose size can be measured. This is a method 
alternative to measuring the frequency of the g c mode in p mode pulsators. 

Several of the SPB and 7 Doradus stars rotate with periods comparable to their oscil- 
lation periods. This, again, calls for models that take rotation into account with a more 
sophisticated approach than perturbation theory; the corresponding work is in progress. 

The 7 Doradus stars are located in a domain where the influence of convection on 
the pulsations is significant; convection is also responsible for the red edge of the 5 Scuti 
instability domain. Modelling with a time dependent convection approach allowed Dupret 
et al. (2005) to reproduce the observed boundaries of these instability strips, and also to 
perform predictions of mode excitation in these stars. The situation for the SPB and 
7 Doradus stars with respect to asteroseismology is therefore the same as for 5 Scuti 
stars: the basic theory is in place, the difficulty remains to find stars permitting the 
extraction of the required information from the observations. 

Beta Cephei stars 

The (3 Cephei stars are massive (~ 9 — 17M Q ) early-B main sequence stars that oscillate 
radially and nonradially in p, g and mixed modes of low radial overtone. This is roughly 
the same type of modes as excited in the 5 Scuti stars, but asteroseismology has been 
more successful for j3 Cephei stars in the recent past due to several reasons. 

The observed frequency spectra are simple enough to provide initial clues for mode 
identification, yet complicated enough to reveal information about the stars' interior struc- 
tures. Photometric and spectroscopic mode identification methods (and combinations of 
both) could be applied successfully to some (5 Cephei stars (e.g., De Ridder et al. 2004), 
additionally aided by the large radial velocity to light amplitude ratios (of the order of 
several hundreds km/s/mag) of their pulsation modes. Sufficient information for unique 
identifications of all modes was obtained; an example was shown earlier. 

Asteroseismic modelling was also eased for /3 Cephei stars, as radial modes have some- 
times been identified; an example was shown earlier. This immediately reduces the pa- 
rameter space in which a seismic model must be sought by one dimension. Due to the 
evolutionary state of j3 Cephei stars (near the centre of the main sequence), a few of the 
observed nonradial modes are of mixed p/g type, which has provided information about 
the convective core size and/or differential interior rotation in a number of stars (e.g., 
Aerts et al. 2003, Pamyatnykh et al. 2004). 

Apart from learning about stellar interiors, asteroseismology of (5 Cephei stars has 
interesting astrophysical implications. Given their high masses, they are progenitors of 
supernovae of Type II which are largely responsible for the chemical enrichment of galax- 
ies. The evolution of massive stars is strongly affected by rotational mixing and angular 
momentum transport (Maeder & Meynet 2000); their internal rotation profile is testable 
by asteroseismology of (3 Cephei stars. 

However, the field has not yet matured to a point where we can claim satisfactory un- 
derstanding of all aspects of the physics governing the interior structures of massive main 
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sequence stars. Several questions may be answered by seismic sounding of (3 Cephei stars: 
how strong is differential interior rotation? How efficient is internal angular momentum 
transport (Townsend 2009)? How strong is convective core overshooting? Can all stars 
between ~ 9 — 17 M & be modelled with the same convective overshooting parameter? 

There are additional questions related to the pulsation physics of j3 Cephei stars that 
need to be addressed. Only the centre of the theoretically predicted j3 Cephei instability 
strip is populated by observed pulsators. Is this an observational shortcoming or a weak- 
ness of theory? What is the upper mass limit of the /3 Cephei stars? Are there post-main 
sequence /3 Cep stars, contrary to theoretical predictions? 

Five different observables strongly depend on the opacities and element mixtures used 
for theoretical modelling of (3 Cephei stars: the radial fundamental to first overtone pe- 
riod ratio, the excited range of pulsation modes, the frequencies of p modes with radial 
overtones larger than two, the dependence of bolometric flux amplitude on the surface 
displacement (Daszyhska-Daszkiewicz & Walczak 2009) and in case of a "hybrid" pulsator 
(see below) the excited range of g modes. Most of these observables are largely indepen- 
dent of each other; modelling of some stars shows that no standard input opacities and 
element mixtures can explain the pulsation spectra in detail. Therefore, the last question 
that may be answered from asteroseismology of j3 Cephei stars is: where must we improve 
the input physics for stellar modelling? 

Some f3 Cephei and SPB stars have emission-line spectra and are thus Be stars. These 
objects rotate rapidly and have circumstellar disks, occasional outbursts etc. They may 
also be studied asteroseismically, but their oscillations are hard to be detected and iden- 
tified, and their rapid rotation requires special theoretical treatment - that is underway. 

Pulsating subdwarf stars 

Three types of pulsating subdwarf star are known: long-period subdwarf B (sdB) stars 
that oscillate in high-overtone g modes (VI093 Herculis stars), short-period sdB stars 
pulsating in low-overtone p and g modes (V361 Hydrae stars), and the only oscillating 
subdwarf O (sdO) star known to date, a low-overtone p mode pulsator. 

Although their g modes would allow the sounding of deep interior regions, their faint- 
ness, long periods and low amplitudes made V1093 Herculis stars escape asteroseismic 
study so far. Theoretical studies of the sdO pulsator have so far been focused on the 
problem of mode excitation. Therefore, the only subdwarf pulsators that have been as- 
teroseismically modelled are among the V361 Hydrae stars. 

This is no easy undertaking as the problem with mode identification and mostly sparse 
(but for a few stars rich and highly variable) frequency spectra again occurs. It is possible 
that these objects have steep interior rotation gradients. Pulsation models must also be 
built upon evolutionary models including the effects of gravitational settling and radiative 
levitation. In practice, modelling is carried out by surveying parameter space in effective 
temperature, gravity, mass and hydrogen mass fraction and by seeking best agreement 
between observed and theoretically predicted oscillation frequencies. Results have been 
obtained on about a dozen of those stars, and a mass distribution consistent with that 
expected from a double star evolutionary scenario has been obtained. Charpinet et al. 
(2009) elaborated on many aspects of asteroseismic modelling of pulsating subdwarf stars. 
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Rapidly oscillating Ap stars 

The rapidly oscillating Ap (or short: roAp) stars are special among the pulsating stars 
because their high-overtone p mode oscillations are predominantly governed by a magnetic 
field that aligns their pulsation axis with the magnetic, and not with the rotation, axis 
of the stars. The magnetic field also distorts the pulsations modes, so they can no longer 
be described with a single spherical harmonic. 

Because of their short periods and very low amplitudes (below 1%) these stars have for 
many years been studied photometrically only; an extensive review was given by Kurtz 
& Martinez (2000). Spectroscopic observations of roAp stars, however, provided a new 
level of insight, particularly into atmospheric physics of these pulsators. 

The reason is that the vertical wavelength of the pulsation modes is shorter than or 
about the same order of the layer thicknesses of the chemical elements in their atmospheres 
that are highly stratified by radiative levitation. Therefore the radial velocity amplitudes 
of the oscillations change with line depth and from one chemical species to the other. 
The chemical elements are also inhomogeneously distributed over the surface, allowing 
three-dimensional studies of the abundances and pulsational geometry (Kurtz 2009). 

Because of the unique possibilities offered by the atmospheric structure of roAp stars, 
spectroscopy is also much more sensitive in detecting oscillations compared to photometry. 
As an outstanding example, Mkrtichian et al. (2008) detected a complete / = — 2 mode 
spectrum for Przybylski's star over 3 — 4 radial orders and performed some initial seismic 
modelling. 

Solar-like oscillators 

The low amplitudes of the stochastically excited oscillations of solar-type stars made 
their observational detection elusive for a long time. In retrospect, the first detection was 
made by Kjeldsen et al. (1995), but confirmed only several years later. Meanwhile, the 
observational accuracy has improved to an extent that detections were made in hundreds 
of stars (mostly giants), and seismic analyses of several were performed. 

The potential of solar-like oscillations for seismic sounding is large. Once detected, the 
pulsation modes are rather easy to be identified because they are high-overtone p modes; 
an example was provided earlier. So are the large and the small frequency separations (if 
detectable), immediately placing main sequence stars on the (asteroseismic) HR diagram. 

Most interesting, as for all pulsators with nearly-asymptotic frequency spectra, are 
irregularities in the latter, caused by features in the stellar interiors. These would for 
instance be the base of their envelope convection zones, or the helium ionization region. 
Houdek (2009) gave an overview of the expected seismic signatures of such features and 
their astrophysical importance. 

The phenomenon of avoided crossing (Fig. 4) also takes place in solar-like oscillators 
once they have reached the subgiant stage, and it makes itself obvious in Echelle Diagrams 
(Bedding et al. 2007). As for other types of pulsator, this would allow an asteroseismic 
determination of the convective core size (for stars sufficiently massive to possess a con- 
vective core). 

The limited lifetimes of the intrinsically damped solar-like oscillations enable inferences 
concerning pulsation mode physics. The observed power spectra at individual mode fre- 
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quencies show a multitude of peaks whose overall shape would correspond to a Lorentzian. 
The half-widths of these Lorentzians yield a determination of the mode damping rates; 
the mode lifetimes are inversely proportional to those. The mode lifetimes are in turn 
dependent on properties of the surface convection zone. 

Theoretical predictions of the amplitudes of solar-like oscillators are important not 
only for understanding their physics, but also for planning observational efforts. After 
years of predictions resulting in amplitudes larger than were observed, the incorporation 
of the mode lifetimes and subsequent computation of mode heights appear to result in a 
scaling law that estimates observed amplitudes well (Chaplin et al. 2009), and seems in 
agreement with measurements up to oscillating giants (Hekker et al. 2009). 

Finally, it is worth to note that as all cool stars possess a convective envelope, it can 
be expected the solar-like oscillations are excited in all of them, up to red supergiants 
(Kiss et al. 2006). 

Hybrid pulsators 

Some of the pulsational instability strips shown in Fig. 1 partly overlap. It is therefore 
logical to suspect that stars that belong to two different classes of pulsating star, having 
two different sets of pulsational mode spectra excited simultaneously, may exist. Indeed, 
a number of those have been discovered. This is good news for asteroseismology as the 
information carried by both types of oscillation can be exploited. The confirmed cases 
have so far always been high-overtone g modes and low-overtone p modes, as evaluations 
of the pulsation constants of the oscillations show. Mixed-mode pulsations by themselves 
are not "hybrid" pulsations because they occupy the same frequency domain as pure 
p modes. 

There are 5 Scuti/7 Doradus stars, (3 Cephei/SPB stars and long/short-period subd- 
warf B pulsators. The main physical difference between the B-type and A/F-type "hybrid" 
pulsators is that in the first group the same driving mechanism excites both types of os- 
cillation, whereas in the 5 Scuti/7 Doradus stars two main driving mechanisms are at 
work. The cooler 5 Scuti stars and all 7 Doradus stars should have thin surface convec- 
tion zones that would support the excitation of solar-like oscillations (e.g., Samadi et al. 
2002). Ongoing searches for such oscillations have so far remained inconclusive, but there 
has been a report of solar- like oscillations in a /3 Cephei star (Belkacem et al. 2009). 

The revision in our knowledge of the solar chemical element mixture (Asplund et al. 
2004) resulted in the prediction of a much larger overlap region between (5 Cephei and SPB 
stars in the HR diagram, and the frequency ranges of the excited long and short period 
modes in "hybrid" B-type pulsators suggest that the heavy-element opacities used for 
stellar model calculations are still too low (Dziembowski & Pamyatnykh 2008, Handler 
et al. 2009). Intriguingly, a similar conclusion has been independently obtained from 
helioseismology (Guzik, Keady, & Kilcrease 2009). 

One of the first hybrid pulsators reported in the literature (HD 209295, Handler et 
al. 2002) turned out to have its g mode oscillations tidally excited by a another star in 
a close eccentric orbit. The changing gravitational influence of the companion gives rise 
to forced oscillations with frequencies that are integer multiples of the orbital frequency. 
The tidal deformation of the pulsating star is similar to a sectoral / = 2 mode, and such 
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modes are therefore most easily excited. There are also several known cases of pulsation 
in ellipsoidal variables. Whereas close binarity represents an additional complication for 
theoretical asteroseismic studies due to the gravitational distortion of the pulsator, some 
other cases of binarity can be used to obtain additional constraints for seismic modelling. 

Pulsation in eclipsing binaries and open clusters 

The fundamental way to determine stellar physical parameters, in particular masses, is 
the analysis of detached eclipsing binary systems whose components can be assumed to 
have evolved as if they were single (Torres et al. 2010). Such constraints are most welcome 
for asteroseismology and therefore it is self-evident that pulsators be sought for in such 
binaries. To date, several dozens of such systems are known. Most of these would be 
5 Scuti pulsators, some are (3 Cephei stars, but the best studied case is a short-period 
sdB pulsator. 

These objects provide another possibility for mode identification (Nather & Robinson 
1974): throughout the eclipse, different parts of the stellar surface become invisible. In 
case of nonradial oscillations, only part of the pulsation modes is seen, and the light 
amplitudes and phases change according to the types of mode. In this way, the oscillation 
mode can be identified, a method known as eclipse mapping or spatial filtration (e.g., 
Reed, Brondel, & Kawaler 2005). 

Another way to gain support for asteroseismic modelling is to study pulsators in stellar 
clusters. Besides the possibility to observe several objects simultaneously (e.g. via CCD 
photometry or multifibre spectroscopy), cluster stars can be assumed to have originated 
from the same interstellar cloud. Therefore, they should be of the same age and chemical 
composition. These parameters can be well determined from the properties of the cluster 
itself, and be imposed as a constraint on the seismic modelling procedure of all pulsating 
cluster members, also known as ensemble asteroseismology. 



A new era m precision 

Most observational results reported in this article so far were based on classical ground 
based observing methods, such as single- and multisite photometry and spectroscopy. 
However, a new era in observational asteroseismology has begun. 

Asteroseismology requires knowledge about as many intrinsic stellar oscillation fre- 
quencies as possible, and these often have low amplitude, calling for measurements with 
the highest accuracy. As this is also a requirement for the search for extrasolar planets, 
synergies between the two fields have emerged. 

Spectroscopically, asteroseismology has benefited from high precision radial velocity 
techniques, such as the iodine cell method, originally invented to find extrasolar planets. 
Measurements of oscillations in distant stars with amplitudes down to 20 cm/s, about one- 
tenth of human walking speed, have become possible in that way. Only about a dozen 
spectrographs worldwide are capable of reaching the required precision. Such observations 
are still therefore expensive in terms of observing time and complexity of data analysis, 
but new observing networks, such as the Stellar Observations Network Group (SONG, 



http://astro.phys.au.dk/SONG/) aim at achieving similar precision on a regular basis. 



32 



Gerald Handler: Asteroseismology 



SIAMOIS (Sismometre Interferential A Mesurer les Oscillations des Interieurs Stellaires, 
Mosser et al. 2009) is expected to work at the same precision and duty cycle with its 
single node placed in Antarctica. 

Concerning photometry, the main problem of ground based observations is scintilla- 
tion, irregular changes in the measured intensity due to anomalous atmospheric refraction 
caused by turbulent mixing of air with different temperatures. The solution to this prob- 
lem is to observe stellar variability from space. Here the synergy with extrasolar planet 
research is that the measurements have sufficient precision to detect the signature of 
transits of planets when they pass in front of their host star. 

Aside from fine-guidance sensor Hubble Space Telescope photometry, the first astero- 
seismic data from space were due to an accident. The main science mission of the WIRE 
(Wide-Field Infrared Explorer) satellite failed due to loss of cooling fluid for the main 
instrument, but the star trackers, of 52 mm aperture, were consequently used for time- 
resolved photometry. An overview of the results can be found in the paper by Bruntt & 
Southworth (2008); WIRE has fully ceased operation in October 2006. 

The first dedicated asteroseismic space mission that was successfully launched is MOST 
(Microvariability and Oscillations of STars, [http://www.astro.ubc.ca /MOST/, Walker et 
al. 2003, Matthews 2006). The spacecraft is in orbit since June 2003 and still continues to 
provide asteroseismic data with its 15-cm telescope as main instrument. One of the most 
interesting results from MOST was the discovery of high-overtone g mode oscillations in 
a supergiant B star (Saio et al. 2006), suggesting an extension of the SPB star instability 
domain to the highest luminosities. 

The CoRoT (Convection, Rotation and Transits, Baglin 2003) mission was successfully 
launched in December 2006 and hosts a 27-cm telescope. It observes two fields of 2 x 3 
degrees on the sky each, one devoted to asteroseismology of a few bright stars, and the 
other searching for planetary transits in many stars, at the same time performing a high- 
precision stellar variability survey. A special volume of Astronomy & Astrophysics (2009) 
reports some of the early CoRoT science. 

Asteroseismically, the detection of solar-like oscillations in almost a thousand giant 
stars is remarkable (Hekker et al. 2009). In addition, one of the most intriguing CoRoT 
results was the detection of large numbers of variability frequencies in two 5 Scuti stars. 
Interpreted as independent pulsation modes, these would supply hundreds of oscillations 
to be asteroseismically modelled, reverting the basic problem for the study of these ob- 
jects: first, there were too few modes available, now there would be too many! However, 
doubts have been raised whether all the frequencies extracted from those data would really 
correspond to normal modes of pulsation, or would rather be a signature of granulation 
(Kallinger & Matthews 2010). 

Originally designed for the search for earth-like planets in the habitable zone, the latest 
addition to the asteroseismic space fleet became Kepler ( ]http: / / www.kepler.nasa.gov] ) , 
launched in March 2009. Asteroseismology can measure the radii of solar-like oscillators, 
among which should be planet-hosting stars, to a relative accuracy of 3%. Thus the radii 
of transiting planets would be known to the same precision. Therefore it was decided to 
devote a small percentage of the observing time of this space telescope with an effective 
95-cm aperture to asteroseismology. Kepler is the most powerful photometry tool for 
asteroseismology to date and observes a ~ 10 x 10 degree field for three years, practically 
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without interruption, providing a time base considerably longer than all other present 
missions. Initial results of the Kepler Asteroseismic Investigation, based on the first 43 
days of science operations, were summarized by Gilliland et al. (2010). 

A possible shortcoming of all these asteroseismic space missions is that they observe 
in only one passband. Therefore all information available for seismic modelling are 
the targets' oscillation frequencies, unless mode identifications are provided by ground- 
based support observations. These are often difficult because the larger space tele- 
scopes mostly observe faint targets. The BRITE-Constellation (BRIght Target Explorer, 
|http: / /www.brite-constellation.at[ ) mission therefore adopts quite a different strategy. 
This mission consists of three pairs of nanosatellites hosting a 3-cm telescope each that 
will observe in at least two different passbands (one per satellite), facilitating mode iden- 
tification with the photometric method. BRITE-Constellation will preferentially observe 
stars brighter than 5th magnitude and has a large 24 x 24 degree field of view. Given the 
brightness of the science targets, mode identification from high-resolution spectroscopy 
can also easily be done. The first pair of satellites are to be launched in early 2011. 

Finally, PLATO (PLAnetary Transits and Oscillations, http://sci.esa.int/plato) is a 
mission designed to provide a full statistical analysis of exoplanetary systems around 
nearby and bright stars. Currently in the definition phase, it will host 28 telescopes of 
10 cm aperture that will observe two 557 square degree fields for 2.5 years each. It is 
intended to observe 100000 stars to a precision of 1 ppm per month and 500000 stars to 
somewhat poorer accuracy to determine stellar and planetary masses to 1%. 

Prospects and problems 

Asteroseismology is a research field evolving so rapidly that some of the results reported 
here will already be outdated when this book appears in print. Given the ~ 30-year 
headstart of helioseismology in comparison, the field now is in its teenage years, but 
matures rapidly. 

The theoretical basis for asteroseismic studies is laid, although far from being perfect. 
Some of the problems that require solution comprise improved treatment of magnetic 
fields, convection, internal flows, and fast rotation. It is still poorly known what causes 
stellar cycles, and what makes certain classes of pulsator select the types of mode they 
oscillate in. Some asteroseismic results point towards a requirement of still higher heavy- 
element opacities, which current calculations do not seem capable of providing. 

Observationally, pulsating white dwarf stars, (3 Cephei stars and solar-like oscillators 
have been studied asteroseismically, and continue to be. There are high hopes that present 
asteroseismic space missions will open the 5 Scuti, SPB, 7 Doradus and V1093 Herculis 
stars for interior structure modelling, and further improve the situation for V361 Hydrae 
stars. Future high-precision radial velocity networks and sites will improve our knowledge 
mostly for solar-like oscillators and roAp stars, with the latter guiding theory of stellar 
pulsation under the influence of rotation and magnetic fields. The Kepler mission will 
provide asteroseismic results for solar-like oscillators en masse, and a large number of 
massive stars are expected to be studied with BRITE-Constellation. The Gaia mission 
:/ /www.rssd.esa.int/gaia[ ) is expected to provide accurate luminosity determinations 
for a vast number of asteroseismic targets, tightly constraining the modelling. 
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It is therefore only appropriate to finish with a quote by Eyer & Mowlavi (2008): now 
is the time to be an asteroseismologist! 
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